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Abstract: — The paper discusses the mathematical model based on the prey-predator system, which diffuses in a near linear 
bounded region and the whole region is divided into two near annular patches with different physical conditions. The model 
has been employed to investigate the population densities of fishes depending on time and position. Finite Element Method 
has been used for the study and computation. The domain is discretized into a finite number of sub domains (elements) and 
variational functional is derived. Graphs are plotted between the radial distance and the population density of species for 
different value of time. 
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I. 1 INTRODUCTION 

The prey-predator relationship is the most well known of all the nutrient relationships in the ocean. [2] Fishes as a 
group provide numerous examples of predators. Thus the herbivorous fishes such as the Sardines and the Anchovy are 
sought after by predators of the pelagic region. These plankton feeding fishes form an important link between the plankton 
and the higher carnivores which have no means of directly utilizing the rich plankton as a source of food. Many of the 
predatory fishes, such as tuna, the barracuda and the salmon, have sharp teeth and they move at high speed in pursuit of their 
prey. [7] Many marine mammals such as the killer whales, porpoises, dolphins, seals, sea lions and walruses are all predators 
with well developed sharp teeth as an adaptation for predatory life. The pelagic fishes of great depths resort to predatory life 
as they have other source of food there. But owing to the scarcity of food available, they are comparatively smaller than the 
predators of the surface region. These fishes have developed many interesting contuvances and food habits. Thus many of 
the deep sea fishes such as the chiasmodus have disproportionately large mouths and enormously distensible stomachs and 
body walls which permit swallowing and digestion of fishes up to three times their own size. The mouth is well armed with 
formidable teeth to prevent the escape of the prey. Some of these fishes are provided with luminescent organs and special 
devices to make the best of the prevailing conditions. Many benthic ferns also are predaceous, living upon each other and 
other bottom animals. [5] Bottom fishes, especially the ray, live on crustaceans, shell fish, worms and coelentrates of the sea 
floor. 



II. MATHEMATICAL MODEL 

In this paper, we discuss the mathematical model based on the prey-predator system, which diffuses in a near 
linear bounded region and the whole region is divided into two patches. The model has been employed to investigate the 
population densities of fishes depending on time and position. [4] Finite Element Method has been used for the study and 
computation. This approach can be extended to the study of population in more patches. 

Variational finite element method is also employed. [1] This method is an extension of Ritz method and is used for more 
complex boundary value problems. In the variational finite element method the domain T is discretized into a finite number 
of sub domains (elements) and variational functional is obtained.[3],[6] The approximate solution for each element is 
expressed in terms of undetermined nodal values of the field variable, as appropriate shape functions (trial functions) or 
interpolating functions. 

Here we have taken a prey-predator model, in which prey-predator populations diffuse between two circular patches in a 
given area. The system of non-linear partial differential equations for the above case is: 
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The total area is divided into two patches. The first patch is assumed to lie along the radius and second 

patch lies along the radius r^ — r — r 2 ■ 
■ th 

Here, for I patch we take 
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=Density of prey population 
P^- =Density of predator population 

=Intrinsic growth rate of prey population 
7*2 • =Intrinsic growth rate of predator population 

=Interspecific interaction coefficient 
}?2i = Interspecific interaction coefficient 

D^ - , D 2i =Diffusion coefficient of prey and predator populations 
Equilibrium points of the equations (1) are 

E u (o,o),E 2i (N*,P*) 
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We take small perturbations U^ - and U 2 - from the non zero equilibrium point i.e. 

n, = n; + u u , p = p* +u 2i ( 2) 
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«1 and U u , U 2i are prey and predator populations. 



Then system of equations (1) becomes 
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Initial conditions are taken as 



«„(r,0) = G 1( (r). 

u 2i (rfi)=G li (r) 

Where G ;i (r) are known functions. 
Interface conditions at T = f\ are assumed to be 
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where and R 2 denote region -I and region-II respectively. 
Boundary conditions associated with the system of equations are 
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III. SOLUTION 

To solve this model, we apply the Finite Element Method. Comparing the system of equations (3) with Euler 
Lagrange's equations, we get 
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The corresponding variational form is 
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Accordingly 
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Substituting values of A ; - , , C ; - , Z) ; - , £" ; - ,0=1.2) and taking Laplace transform of (15)and (16), we get 
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Here 1 j ,U 2 2\ are Laplace transforms of Wjij , ^221 res P ect i vei y and U^^q,U^^ 2 ,U 22 q,U 222 are independent 
of t and Mi i j (0) and W221 (0) are initial values. Solving (17 and (18), we get, 
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Taking inverse Laplace transform of 1 j and U 22 ^ ,we obtain the expressions for and U 22 \ .Substituting these 

values in (13), we get values of and U 2i (i=l,2).Using these values we can get the values of and P- (i=l,2). 

IV. NUMERICAL COMPUTATION 

We make use of the following values of parameters and constants. 
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r 12 = .04 , b n = .5 , b n = A , K n = 100 , K n =125 , d n = .9 , d 12 = .8 , r 21 = .025 , 
r 22 =.035 ,fc 21 = .5,6 22 =.4,^ 21 =110 ,K 22 =130 ,d 21 =.8, d 22 =.9,r 0 =10 =20 , 
r 2 = 30 ,m 1110 = .5 w 2210 = .5 ,w 1112 = .9,w 222 = .9,m 110 = .2,w 220 = .2. 

V. CONCLUSION 

Graphs are plotted between radius r and population density of species for different values of time. Graphs show 
that density of prey increases with radius while density of predator decreases with radius at any time in the first patch. In the 
second patch density of prey decreases with radius while density of predator increases. Graphs also show that density of 
prey increases with time and density of predator decreases with time in the first patch. In the second patch density of prey 
decreases with time while density of predator increases with time. 
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